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Abstract 

At upcoming medium baseline reactor neutrino experiments the spallation ®Li background 
will be somewhat larger than the inverse ft decay reactor neutrino signal. We use new 
FLUKA simulations of spallation backgrounds to optimize a class of veto strategies and 
hnd that surprisingly the optimal veto for the mass hierarchy determination has a rejection 
efficiency below 90%. The unrejected background has only a modest effect on the physics 
goals. For example for the hierarchy determination falls by 1.4 to 3 points depending on 
the muon tracking ability. The optimal veto strategy is essentially insensitive to the tracking 
ability, consisting of 2 meter radius, 1.1 second cylindrical vetoes of well tracked muons with 
showering energies above 3 to 4 GeV and 0.7 second full detector vetoes for poorly tracked 
muons above 15 to 18 GeV. On the other hand, as the uncertainty in 612 will be dominated 
by the uncertainty in the reactor neutrino spectrum and not statistical fluctuations, the 
optimal rejection efficiency for the measurement of 612 is 93% in the case of perfect tracking. 
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1 Introduction 


The next generation of medium baseline reactor neutrino experiments [1, 2] will deliver by far 
the most precise measurement yet of Qyi [3] and also will be sensitive to the neutrino mass 
hierarchy [4] . Cosmogenic muons interacting with in their organic liquid scintillator 
detectors will create spallation products ®Li and ®He whose decays yield the same double 
coincidence used to detect reactor neutrinos [5]. The background will be larger than the 
signal [6]. The effectiveness of a veto strategy depends critically on the ability to track 
these muons, so that only events near the track need be vetoed. In this note we present the 
optimal strategy as a function of the tracking efficiency. Our main result is that the optimal 
strategy leads to the acceptance of a much larger fraction of the background than had been 
found in previous studies but that this large background has a relatively minor effect on the 
main science goals of these experiments. 

On a smaller scale, this issue has already been confronted by KamLAND [7]. As the 
detector is deeper and smaller, the muon rate is only 0.2 Hz, appreciably slower than the 
257 ms decay time of ®Li. Of these muons, the most dangerous are the showering muons. 
These are defined to be muons which deposit at least 3 GeV in the detector in addition 
to the energy deposited by ionization. While they occur at a rate of only 0.03 Hz, they 
are responsible for the majority of the ®Li background [8]. As a result, KamLAND chose a 
full detector veto of 2 seconds after all showering muons and also after muons whose track 
could not be reconstructed. For nonshowering muons whose track could be reconstructed, 
a 3 meter radius cylinder about the track was vetoed for 2 seconds. The combination of 
these vetoes left a small number of background events, estimated for example as 14 [9] of 
the original 10^ [7]. Overall, this strategy led to a rejection efficiency of about 99% while 
removing only about 15% of the signal. 

Can this same strategy be repeated with upcoming experiments like JUNO and RENO 
50? JUNO will have a muon rate of about 5 Hz and a showering muon rate of about 1 
Hz, and so a 2 second veto after each showering muon would eliminate essentially all of 
the events [6]. However we will show in this paper that a much less efficient veto can 
give an acceptable dead time and result in a relatively modest effect on the main science 
goals of the experiment, the determination of the neutrino mass hierarchy and the most 
accurate ever measurement of Qyi- The reason for this is that the shape and normalization 
of the ®Li decay spectrum can be determined fairly well, and the ®He decays yield a small 
background. Therefore the background can be subtracted, leaving only a small contribution 
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to the statistical fluctuations. In this note we will determine the optimal veto strategy and 
its impact on the sensitivity to the hierarchy and the measurement of Qyi- 

This procedure is somewhat complicated by muon bundles. At both KamLAND and 
JUNO, nearly half of the muons that hit the detector are contained in bundles of multiple 
muons created in the same cosmic ray collision. However, as the KamLAND detector is 
much smaller than the characteristic muon separation, in less than 5% of events do multiple 
muons from the same bundle strike the detector. On the other hand, in the case of JUNO, 
about 10% of muon events will consist of multiple muons striking the detector. The tracking 
of these bundle events is extremely challenging, and so the ability to veto only the muon 
track, and not the full detector, will depend on whether JUNO is capable of tracking multiple 
muons. In this study we will examine the dependence of the sensitivity to the hierarchy and 
the precision with which Qyi can be measured upon the tracking ability of the detector. 

2 The ®Li Background 

51% of ®Li and 16% of ®He decays yield a neutron together with an electron, mimicking 
the double coincidence signal of inverse (3 decay (IBD). At JUNO these double coincidence 
background rates will be 8.5 x 10“^ Hz and 3.0 x 10“® Hz respectively. The two decay spectra 
are so similar that so far KamLAND has not been able to conclusively observe ®He beyond 
the 1 a level. Therefore, to a reasonable approximation the ®He background may be ignored. 

2.1 The Distribution of ^Li 

At KamLAND a full detector or cylindrical veto was applied to all cosmogenic muons. 
However those muons whose showering energy, defined to be the energy that they deposit in 
the scintillator minus 1.43 MeV per cm of their track^, is sufficiently small are exceedingly 
unlikely to generate any ®Li. In fact, we will see below that in the case of JUNO, there are so 
many muons that the optimal strategy is not to apply any veto to muons whose showering 
energy is less than a certain value, which is at least about 3 GeV in each case. The average 
muon track length is about 24 meters. To determine the volume-weighted dead time, we will 
need the average track length of only those muons whose showering energies exceed a certain 
minimum Ef. We have found that this average track length is 28 meters whenever Et > 1.5 
GeV, as it for the veto strategies that we will find to be optimal. Therefore, while we will 

^1.43 MeV/cm is the ionization energy deposited in electrons with less than 100 keV. Those with more 
energy are considered to be 6 rays and little would be gained by subtracting them as well from the showering 
energy [6]. 
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Figure 1: The radial distribution (horizontal) of ®Li about a muon track, normalized to unity 
for each muon showering energy (vertical) 


use the average track length d{Et) as determined by our FLUKA simulations, in practice for 
the optimal veto strategies this will always be equal to 28 meters. 

The distribution of ®Li events about this track then depends on the muon’s showering 
energy and also the distance from the muon track. As can be seen in Fig. 1, this distribution 
to a very good approximation factorizes, in the sense that the radial distribution of ®Li is 
essentially independent of the showering energy. This greatly simplihes our analysis, as it 
means that the ®Li yield is entirely characterized by just two functions, the radial, shown in 
Fig. 2 and showering energy distribution which was already provided in Ref. [6]. 

In the next sections we will attempt to determine the optimal veto strategy. To do this, 
we will need the fraction of muons with a showering energy above a threshold and the 
fraction of ®Li produced by those muons. These are displayed respectively in Figs. 3 and 4. 
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Figure 2: The radial distribution of ®Li about a muon track 



Figure 3: The fraction of muons with a showering energy above the threshold Et GeV 



Figure 4: The fraction of ®Li generated by muons with a showering energy above the threshold 
Et GeV 
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3 Live Time vs Rejection Efficiency 


3.1 The veto strategy 

In practice the optimal veto strategy may be quite complicated, and any strategy should be 
rehned once data taking begins. Our goal in this paper is to understand the effect of the 
spallation isotope backgrounds on the science goals of medium baseline reactor experiments 
and the implications of these backgrounds for the muon tracking requirements. As a result, 
it is necessary to produce a reasonably effective but simple veto strategy, which can be used 
in the analysis. 

The strategies that we will consider consist of two cuts 

1) When the energy E > Ej and it is known that the muon (bundle) was not tracked 
successfully, veto the full detector for a time tf. 

2) When the E > Ec and it is believed that the muon was tracked successfully, veto 
a cylinder of radius Tc about each muon track for a time tf.. 

We will optimize Ef, Ec, Tc, tf and tc separately for the neutrino mass hierarchy deter¬ 
mination and for the measurement of 6 * 12 . We will first assume that each muon is tracked 
perfectly. This assumption will be relaxed in Sec. 5. 

One might think that the simultaneous optimization of these 5 variables would require a 
slow, global search through a 5-dimensional space, rerunning our simulations at each step. 
In fact the simulations only depend on two parameters, the background rejection efficiency 
and the live time. Therefore we will use the following strategy. First, in this section, using 
the results of Sec. 2.1 we will find the strategy which maximizes the live time for each 
possible background rejection efficiency. This maximum is independent of the science goal 
considered, and so this section will apply to both the hierarchy search and to 9 i 2 , and can 
in the future be applied to other potential science goals such as the study of geoneutrinos 
should the reactors all be turned off in the future. 

Next, in Secs. 4 and 5 simulations of the detector JUNO with respectively perfect and 
imperfect tracking are considered with each possible background rejection efficiency and live 
time. As the optimal live time for each detector efficiency will already be determined, for 
each science goal we will only need to optimize a single variable, the background rejection 
efficiency. This will yield the optimal cuts as well as the effect of the background upon the 
hierarchy determination and the measurement of 612 . 
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Let m{E) and 1{E) respectively be the fraction of single mnons and the fraction of ®Li 
isotopes created by single mnons with showering energy greater than E. These two fnnctions, 
in the case of the detector JUNO are plotted in Figs. 3 and 4. Here we have not considered 
mnon bnndle events. While only 10% of mnon events at JUNO will be mnltimnon events, 
these will acconnt for abont a third of all mnons and so abont a third of spallation isotope 
prodnction. 

In the case of mnon bnndles, the spallation isotope prodnction rate is jnst the snm of the 
rates for the individnal mnons. The energy distribntion for mnon bnndles is similar to that 
of single mnons and so an extension of onr analysis to bnndle events wonld be straightfor¬ 
ward. One need only modify the showering energy distribntion m{E) and spallation isotope 
distribntion 1{E) fnnctions to be the weighted snm of the distribntions for all fc-mnon bnn¬ 
dles, where the total deposited energy is the total of that of the k mnons. This reqnires no 
new simnlations, as the fc-bnndle distribntion is essentially jnst the snm of k single mnon 
distribntions. One only needs to know the probability of mnltimnon events, which in the 
case fc = 2 is calcnlated in Ref. [6]. 

The inclusion of bundles would reduce the rejection efficiency as the scintillator detectors 
cannot determine how much energy was deposited by each muon, and so one would need 
to rely upon a total energy threshold. As a result, the rejection efficiencies obtained in this 
study are optimistic estimates. However, it will be clear below that even a 10% reduction in 
the rejection efficiency would not have an enormous impact on the science goals. 

3.2 Full detector veto 

A full detector veto of time tf after the passage of each muon with a showering energy of at 
least Ef will yield a rejection efficiency of 

£,{E,,t,) = 1{E,){1 - (3.1) 

where ~ 255 ms is the decay time for ®Li. Only {3 decays of ®Li to excited states of ®Be 
eventually yield free neutrons and so a false double coincidence, however the decay of ®Be 
into two a particles and a neutron is much faster than the (3 decay of ®Li and the (3 decay 
time of ®Li is independent of whether the ®Be is in an excited state, and so for the decays 
producing neutrons is equal to the overall tf for ®Li. 

For a given rejection efficiency Sf, one can invert (3.1) to obtain Ej as a function of the 
veto time tf 

B/(0.</) = r>(r-L). (3.2) 
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Figure 5: The optimal full detector veto time tf{Sf) yielding an ®Li rejection efficiency Sf 
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Figure 6: The optimal full detector showering energy threshold Ef{Sf) and fractional live 
time ti{Sf) yielding an ®Li rejection efficiency £f 


The resulting fractional live time ti is determined from the product of the veto rate m{Ef) 
and the veto time tf 

tf) = exp {-tf m{Ef{Sf, tf))). (3.3) 

The veto may then be optimized, for each Sf,hj simply setting 

= 0. (3.4) 

dtf 

The solution to this equation is the optimal tj for each rejection efficiency Sf, which we will 
call tf{Sf). The optimal veto time tf{Sf) for each efficiency Sf is shown in Fig. 5. Note that 
the optimal veto time is about 1 second, much less than the 2 seconds at KamLAND. 

Once tf{Sf) is found numerically, Ef{Sf) is found using Eq. (3.2) and the fractional 

live time ti{Sf) using (3.3). These are shown in Fig. 6. Note that even a modest rejection 

efficiency of 90% requires a showering energy threshold of about 4 GeV and leads to a live 
time of less than one half of the total runtime. In a 20 kton detector, full detector vetoes 
are quite costly. As a result, we will assume that they are only used when the muon was not 
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tracked successfully, a case which will be considered in Sec. 5. In Sec. 4 we will assume that 
all muons are well tracked and so only cylindrical vetoes will be considered. 


3.3 Cylindrical veto 


The efficiency of a cylindrical veto depends not only on the functions m{E) and 1{E) char¬ 
acterizing the muon distribution and ®Li yields as functions of energy E, but also upon the 
®Li radial distribution about the muon track f{R). Here f{R) is the ®Li rate per unit radius 
at a radial distance R from the track, as illustrated in Fig. 2. A radius Tc cylindrical veto 
of time tc after the passage of each muon with a showering energy of at least Ec will yield a 
rejection efficiency of 

£AE,,tc, To) = i(i5o)(l - (3.5) 

Jo f{R)dR 

For simplicity we will dehne the integral 


E{rc) 


/;• f{R)dR 

r 


which is equal to the normalized ®Li event rate at radii below Vc- 


(3.6) 


Note that the expression (3.5) for the rejection efficiency Tc factorizes into a part depen¬ 
dent upon Ec and a part x dependent upon tc and Vc 


Sc{Ec,tc,rc) = l{Ec)x{tc,rc), 


x(tc, fc) = (1 - e F(rJ. 


(3.7) 


Therefore we may separate the optimization problem into steps. First we may £x x. This 
gives a relation between tc and Tc 

tcix, rc) = -td In (^1 - ■ (3-8) 

Now instead of the three independent variables {Ec, tc, Tc) there are three independent vari¬ 
ables {Ec,rc,x). However the rejection efficiency Sc is entirely determined by Ec and x via 
Eq. (3.7) and so Vc can be varied, hxing Ec and x, so as to minimize the total dead time. 
Here the volume-weighted dead time that is relevant is the veto time multiplied by the veto 
volume corresponding to a live time of 

ti{Ec, Tc, x) = exp rc)m{Ec)^ (3.9) 

where V is the detector volume and, following the dehnition in Sec. 2.1, d{Ec) is the average 
tracklength of a muon with showering energy Ec- The detector efficiency will be proportional 
to ti. 






Fc (meters) (seconds) 




Figure 7: The optimal cylindrical veto radius Tc and time tc for a fixed value of x, defined 
in Eq. (3.7) 


We want to choose so as to maximize ti. However only enters in the combination 
r‘^tc{x,rc) and so we only need minimize r%ix,rc) for each value of x 


dr^tcix, 



(3.10) 


Using Eq. (3.8) this minimization is easily performed, yielding the optimal value of rc{x) 
and therefore tc{x) for each value of x. The optimal values are shown in Fig. 7. Similarly 
to the full detector vetoes, the optimal veto time is about 1 second, only half as long as 
the KamLAND vetoes. Similarly the optimal veto radius is always less than the 3 meter 
KamLAND veto. Due to KamLAND’s greater depth and smaller size, the muon rate is much 
lower which allows KamLAND to use longer vetoes with a small effect on the measurement of 
the reactor mixing parameters. While an analysis such as ours applied to KamLAND would 
nonetheless find that the optimal veto would be longer than at JUNO but shorter than 2 
seconds, nonetheless as a result of the low muon rate KamLAND does not need a precise 
optimization of the veto time. On the other hand, at JUNO the conflicting constraints 
imposed by the large background, which requires a long veto, and the statistical fluctuations 
in the hierarchy-sensitive region of the spectrum around 3 MeV, which require a short veto, 
leave a small window of acceptable veto strategies for the hierarchy determination at JUNO. 

We have reduced the number of independent parameters from 3 to 2, only E^. and x need 
to be optimized. The problem is now identical to that of the full detector veto treated above 
but now X plays the role that was played by tf, or strictly speaking by 1 — Fixing 

the rejection efficiency Ec-, Eq. (3.7) gives the showering energy threshold E^ 


Ec{£c,x) = I ^ 



(3.11) 
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Figure 8: The optimal value of x, defined in Eq. (3.7), as a function of the efficiency Ec of a 
cylindrical veto 


Substituting this value of Ec{£c,x) into (3.9) one obtains ti{Ec,x). Thus, as in the case of a 
full detector veto (3.4), for each rejection efficiency Ec one determines the optimal x{Ec) by 
demanding that it extremizes ti{Ec,x) with Ec held fixed. This optimum is shown in Fig. 8. 

Once x{Ec) has been determined, Eq. (3.11) yields the threshold showering energy Ec for 
the cylindrical veto, shown in Fig. 9. For an efficiency Ec near 99% this threshold should be 
quite low, near 1 GeV. However a more modest efficiency of 90% requires a threshold of 3 
GeV, equal to the threshold used by KamLAND for its full detector veto. 

Gombining Fig. 8 for a; as a function of Ec with Fig. 7 which provides Tc and G as a function 
of X one finds the cylindrical veto parameters Vc and tc as a function of the efficiency £c- 
These are shown in Fig. 10. Finally, as the fractional live time is now quite close to unity, 
we plot the fractional dead time 1 — in Fig. 11. This is one of the main results of our 
study. Note that a rejection efficiency of 90% can be obtained while losing only 3% of the 
IBD signal, while a 95% efficiency costs 10% of the signal. However the commonly quoted 
goal of 99% background rejection would require a loss of 35% of the IBD events, and this 
under the assumption of perfect tracking which is quite unlikely for 99% of muons. 

4 Optimal Veto for Each Science Goal 

In Sec. 3 for each kind of veto, full detector or cylindrical, we obtained the maximum 
fractional live time ti possible for each ®Li rejection efficiency Ef or Ec and we found the 
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Figure 9: The optimal showering energy threshold as a function of the efficiency 8 ^ of a 
cylindrical veto 
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Figure 11: The fractional volume-weighted dead time 1 — t;, as a function of the efficiency 
Sc of a cylindrical veto 


corresponding veto parameters Ef, Ec, tf and tc which are functions of the rejection 
efficiencies. 

In this section we will assume that the muon tracking is perfect, leaving the general case 
to Sec. 5. Therefore we will only be interested in cylindrical vetoes. The cylindrical veto is 
then fully described by a single parameter, the efficiency Sc- The showering energy threshold 
Eci.Sc), veto cylinder radius rdSc) and veto duration tc(i’c) may be read from Figs. 9 and 10 
while the corresponding fractional dead time was reported in Fig. 11. 

But what is the optimal value of Sc^ This depends on the science goal considered. 
For each science goal, we have calculated a statistic characterizing the goal. This statistic 
depends on the total live time of the experiment and on the total background during that live 
time. We have assumed that the experiment runs for 6 years, and so the live time will just 
be Qti years. As ti is known as a function of Sc (Fig. 11), one may then obtain the statistic 
as a function of Sc alone and then choose the Sc which optimizes the statistic, optimizing 
the veto strategy for that science goal. As the vetoes are not hardwired, there is no problem 
using a different veto strategy for each science goal. 

4.1 Neutrino mass hierarchy 

In this subsection we will optimize the veto strategy for the determination of the neutrino 
mass hierarchy at JUNO. For this purpose, we have computed the statistic which is 
equal to the difference between the value of a best fit of a given spectrum to the normal 
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Figure 12: The Ax^ with which the neutrino mass hierarchy may be determined by JUNO in 
6 years with a fraction ti of live time and ®Li rejection efficiency Sc- For the curves from top 
to bottom the fractions ti are 1.0, 0.9, 0.8, 0.7, 0.6 and 0.5. It is assumed that the nonlinear 
energy response of the detector is known perfectly. 


and inverted hierarchies 

(normal) — y^ (inverted). (4.1) 

In Fig. 12 we display Ay^ computed for an Asimov data set generated using the inverted 
hierarchy, using 6 years of run time, which throughout this paper we assume corresponds to 
1.35 X 10® inverse (3 decay events. 

Note that Ay^ is not that of Wilks’ theorem, which is dehned to be the difference between 
a best y^ £t to a null hypothesis and that to a best £t of a continuous family of parameters 
of which the null hypothesis is one choice. Indeed, the choice of hierarchies is not even 
continuous. For example, the quantity Ay^ dehned in Eq. (4.1) is Gaussian-distributed 
[10, 11] whereas that in Wilks’ theorem is distributed following a y^ distribution. Therefore 
this Ay^ is not the square of the number of a of conhdence in a null hypothesis. However, 
in a Bayesian framework in which one assigns a symmetric prior to the two hierarchies, this 
Ay^ is simply related to the sensitivity s to the hierarchy via the relation [10, 11] 


1 

l + e-^xV2- 


(4.2) 


If this sensitivity is converted into a number of a using the error function, one hnds a 
reduction of about 0.5cr as compared with the naive square root formula. 

The detector energy resolution, number of events, etc. are as in Ref. [12] and conform 
to the experimental goals. However we assume that the nonlinear energy response of the 
detector is perfectly known. In fact the experimental goal is a 1% precision on the nonlinear 
energy response, which has a large effect on Ay^ as has been foreseen in Ref. [13] and 
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Figure 13: The Ax^ with which the neutrino mass hierarchy may be determined by JUNO 
in 6 years as a function of the rejection efficiency Sc- The optimal rejection efficiency is 90%. 


demonstrated in Ref. [14]. However the precise effect on depends heavily on additional 
assumptions, such as the shape of the nonlinearity. It may be possibile to estimate the 
nonlinearity by using the known general features of the spectrum [1] , however such a strategy 
may be ineffective or even problematic if the reactor neutrino spectrum has unexpected 
features, such as the 5 MeV bump hrst reported by RENO [15] and later conhrmed by 
Double Chooz [16] and Daya Bay [17]. Therefore the inclusion of the unknown nonlinear 
energy response in this study would introduce uncertainties which would obfuscate the effect 
of the spallation backgrounds. 

The optimization of the rejection efficiency Sc is now quite easy. One need only substitute 
Fig. 11, which gives ti as a function of Sc, into Fig. 12 which gives Ay^ as a function of ti 
and Sc to obtain Ay^ as a function of Sc alone, shown in Fig. 13. One can see that, in 
contrast with the common claim that the optimal rejection efficiency at JUNO will be at 
least 99%, the optimal rejection efficiency is in fact only 90% when the tracking is perfect. 
Using Figs. 8, 9, 10 and 11 with Sc = 0.9, one may obtain the corresponding optimal veto 
parameters, which are summarized in Table 1. 

Despite the relatively large contamination, in the case of perfect tracking the ®Li back¬ 
ground only reduces Ay^ from 13.8 to 12.4. This is our main result. 

If Taishan reactors 3 and 4 are not built, then the reactor neutrino signal at JUNO will 
be reduced by about 26%. This can be incorporated into the present analysis with a 26% 
reduction of the live time ti in Fig. 12, resulting in a 20-25% reduction in Ay^. 
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4.2 Uncertainty in the ®Li spectrum 

The small reduction in is a result of the fact that the shape of the ®Li decay background 
is known quite well, and so can be reliably subtracted. What is the shape? First, recall that 
we are only interested in decays producing neutrons 

®Li —)■ ®Be* + e“ + W —)■ n + 2 q! + e“ + W (4.3) 

where ®Be* is an excited state. The W is not observed. The spectrum of the e~ is known quite 
precisely from Fermi’s theory of /3 decay. It extends up to about 11 MeV and overlaps with 
the entire reactor neutrino prompt energy spectrum. Indeed the scintillator detector cannot 
distinguish these e~ from the e'^ produced by the inverse /3 capture of reactor neutrinos. 
The neutron and a particles are produced essentially immediately after the (3 decay. 

However, as 90% of neutron-emitting decays begin with a (3 decay to the low energy 
(less than 3 MeV) excited states of ®Be [18], little energy is available for the neutron and a 
particles. As a result, these particles have kinetic energies in general well below 2 MeV, often 
below 1 MeV. In this regime their scintillator quenching factors are large. The a particle 
quenching factor is more than 10 times that of the electron, reducing its contribution to 
the scintillation light by a factor of more than 10 so that its contribution to the observed 
energy negligible. The neutron’s quenching factor is more than 5 times that of the electron, 
so its contribution is small but not negligible. The neutron’s energy begins to arrive in the 
scintillator as it slows. The neutrons slows as it collides with particles in the scintillator. 
The time between such collisions is about 5 ns and about half of the collisions, those with 
free protons, reduce the neutron’s energy by about one half. Therefore the neutron loses 
most of its energy in about 10 ns. 

While one might think that the PMTs could then distinguish the neutron’s energy from 
that of the electrons using the 10 ns time difference and so veto the ®Li events, this is not the 
case. The main problem is the time-lag before the scintillator scintillates, which ranges from 
10s to 100s of ns. Another problem is that JUNO is about the same size as the mean free 
path of the ultraviolet photons created by the wave shifters. As a result, a significant portion 
of the photons created by the electron bounce inside of the scintillator before arriving at the 
PMTs, and so arrive 10s of nanoseconds late. These late photons from the electron cannot be 
distinguished from the photons from the neutron slowing, and so one only observes the total 
energy. However, the neutron energy was already 5-10 times smaller than the electron energy 
and the quenching factor reduces its contribution to the visible energy by an additional factor 
of 5, thus the neutron contributes only about 5% of the visible energy. Now the neutron 
spectrum, although it has been measured by many experiments [19, 20, 21], is not known 
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as precisely as the electron spectrum. Nonetheless, the fact that it contributes only 5% of 
the visible energy means that the uncertainty in the neutron spectrum yields a very small 
uncertainty in the ®Li spectrum. This is why the ®Li spectrum is so well known, which allows 
the presence of a large ®Li contaminant in the IBD event sample to have such a small effect 
on the sensitivity to the hierarchy. JUNO itself will provide a very precise measurement of 
the quenching factor weighted spectrum, actually of a weighted combination of the ®Li and 
®He spectra, which will only decrease this uncertainty. 

In fact, the biggest uncertainty in the ®Li spectrum may arise from the veto efficiency 
Tc, which depends upon the the precision of the reconstruction of the muon track and the 
radial distribution f{R) of the ®Li which is not well known, although JUNO will measure 
it better than it has ever been measured before. Unlike the previous uncertainty, this does 
not affect the shape of the observed energy spectrum, but only the overall normalization. 
The precision of these estimates are quite difficult to forecast. Therefore in this note we 
simply assume that the shape of the ®Li spectrum is known perfectly and we include a 1% 
uncertainty in its overall normalization. 

4.3 Mixing angle Ou 

In this subsection we apply the same analysis to the determination of 6 ^ as that used for 
the determination of the hierarchy in Subsec. 4.1. We will choose the veto strategy to 
minimize cr(sin^(26'i2)), the precision with which sin^(26'i2) can be determined. There are 
various calculations of this uncertainty in the literature. The earliest [22, 23] do not consider 
the uncertainty in the shape of the reactor neutrino spectrum. Later calculations, such as 
Ref. [24] , do consider this uncertainty but use the theoretical lower bounds on the uncertainty 
given in Ref. [25] . However, it is now known that the lower bound on this uncertainty is far 
from saturated, indeed it would imply that the 5 MeV bump [15, 17] is excluded at 4(T. 

This motivates us to use the more recent determination of a{siv?{26i2)) from Ref. [26], 
shown in Fig. 14, which uses the observational uncertainty on the reactor neutrino flux from 
Daya Bay [17]. In addition, this is the only study which yields (j(sin^(26'i2)) as a function of 
the dead time and the ®Li veto efficiency, which is necessary to optimize the vetoes. 

The most obvious distinction between Fig. 14 and the corresponding hgure for the mass 
hierarchy. Fig. 12, is that in the former the experimental performance hardly depends upon 
the rejection efficiency and the dead time in the range plotted. This is because the precision 
of a determination of di 2 is dominated by systematic errors, the uncertainty in the reactor 
spectrum, unlike the determination of the hierarchy. Of course, if we included the uncer¬ 
tainty in the detector’s energy response then then the systematic errors would play a more 
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cr(sin^(20i2)) 



&c (Rej. Eff.) 


Figure 14: The Icr uncertainty with which sin^(26*i2) can be measured at JUNO in 6 years 
with a fraction ti of live time and ®Li rejection efficiency Sc, from Ref. [26]. For the curves 
from bottom to top the fractions ti are 1.0, 0.9, 0.8, 0.7, 0.6 and 0.5. The main contribution 
to the uncertainty arises from the uncertainty in the reactor spectrum measured by Daya 
Bay [17]. 


important role in the determination of the hierarchy. However, so long as this response can 
be understood well enough to achieve a reasonable determination of the hierarchy, it will 
play a subdominant role. 

As in the case of the hierarchy, we now substitute Fig. 11, which gives ti as a function of 
Sc, into Fig. 14 which gives cr(sin^(26'i2)) as a function of ti and Sc to obtain a{sm‘^{29i2)) as 
a function of Sc alone, shown in Fig. 15. Note that the optimal rejection efficiency is now 
93%, as compared with 90% for the mass hierarchy and the precision cr(sin^(26'i2)) is 0.0038 
or 0.45%, as compared with 0.44% with no background at all. As summarized in Table 1 this 
means that the optimal veto strategy for the determination of 9 i2 should have 6% dead time, 
as compared with only 3% for the neutrino mass hierarchy. The reason for this is that the 
error in 9 i2 is dominated by the systematic error in the reactor spectrum, not by statistical 
fluctuations like the measurement of the hierarchy. Therefore a larger dead time in the 9 i2 
analysis has less effect on the precision of the measurement than it would in a determination 
of the hierarchy. 

Physically this is a consequence of the fact that 9 i2 is determined by a broad feature in 
the reactor neutrino spectrum, essentially the difference between the neutrino flux at inter¬ 
mediate energies at the solar oscillation maximum and that at higher and lower energies. 
Being a broad feature it is insensitive to the energy resolution of the detector and more 
importantly it involves essentially the whole sample of neutrinos. Thus the statistical un¬ 
certainty on 9 i2 is extremely small. The determination of the hierarchy is quite different. 
It comes from comparing AM^g [27] measured at intermediate and high energies with the 
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o-(sin^( 26 'i 2 )) 
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Figure 15: The la uncertainty with which sin^(20i2) can be measured by JUNO in 6 years as 
a function of the rejection efficiency £c- The optimal rejection efficiency is 93%, higher than 
for the hierarchy because the uncertainty is dominated by a systematic error, the uncertainty 
in the reactor neutrino spectrum, which is taken to be as in Fig. 14. 


Science Goal 

Statistic 

Tc 

ti 

Ec 

X 

tc 

Tc 

Mass hierarchy 

= 12.4 

90% 

97% 

2.9 GeV 

95% 

1.1 sec 

2.3 m 

Measure 612 

<7(sin"(2ei2)) 

93% 

94% 

2.4 GeV 

96% 

1.2 sec 

2.5 m 

: 2 yj.^o/0 

sm ( 2012 ) 


Table 1: The optimal cylindrical veto parameters for the determination of the mass hierarchy 
and di 2 assuming perfect tracking 


position of the 15th, 16th and 17th peaks at low energies [28, 23]. These low energy peaks 
contain few events, especially when one considers the effect of the hnite energy resolution 
and interference from neutrinos coming from reactors at distinct baslines [28] at JUNO. 
Therefore statistical uncertainty is very important in the determination of the hierarchy, 
which drives the preference for a veto strategy with little dead time. 

5 Imperfect Tracking 

In this section we will consider the effects of imperfect tracking of muons. We will make 
the poor approximation that the probability of tracking a muon is independent of its shower¬ 
ing energy. Up to this point the generalization to muon bundles has been straightforward, it 
only required modifying the functions m{E) and 1{E) by including the weighted sum of bun¬ 
dles. However if one introduces bundles, then one either needs to introduce new parameters 
describing how well the bundles are tracked, or else make the further crude approximation 
that the tracking of each muon is independent of whether it arrives in a single or multimuon 
event. In the latter case, the treatment below remains applicable, one need only remem- 
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Variable 

Tracking 

1 — a — (3 

Muon well tracked 

a 

Muon poorly tracked but believed to be well tracked 


Muon poorly tracked and known to be poorly tracked 


Table 2: The tracking of muons is characterized by the parameters a and (3, assumed to be 
independent of the showering energy 

ber that a and f3, which will be dehned momentarily, are fractions or muons, not of muon 
bundles. 

With these caveats, the tracking is described by just two parameters, a and /?, which, as 
is summarized in Table 2 are the fraction of muons for which the tracking fails but is believed 
to succeed and the fraction of muons for which the tracking fails but the experimenter realizes 
that the tracking has failed. Our veto strategy then demands that full detector vetoes be 
applied to a fraction (3 of all muons. In the case of the a muons a cylindrical veto is applied, 
but as it is quite unlikely that this cylinder coincides accidentally with the muon track, we 
will simply assume that the veto does not increase the rejection efficiency. The fraction of 
muons for which tracking is successful is 1 — a — (3. 

For simplicity we will make the crude approximation that the tracks of well-tracked muons 
are determined precisely. In fact, one expects a precision of only 1-2 meters. However, this 
precision depends on the conhguration of the top veto, which is not yet fixed, and so cannot 
be estimated reliably at this time. The result of the uncertainty in the muon track position 
would be to increase the optimal Tc by 1 meter or more and indeed it may be worthwhile 
to consider a cylinder which gets wider along the muon’s trajectory. As a consequence, the 
total dead time will be higher than that estimated below and the rejection efficiency lower. 

The resulting fractional live time is 

ti{Ec,rc,x) = exp |^-(1 _/j) tcm{E^) - /3tfm{Ef)^ . (5.1) 

The veto strategy now depends on two parameters, the rejection efficiencies Sf and Sc for the 
full detector and the cylindrical vetoes. On the other hand, the science goals only depend 
upon ti and the total rejection efficiency 

S = {l-a- (3)Sc + ^Sf. (5.2) 

Therefore the optimization of the veto strategy will consist of two steps. First, for each S, 
we will hnd Sc and Sf so as to maximize ti. Next, in the following two subsections, we will 
choose S for each science goal so that the corresponding statistic is optimized. 
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Figure 16: The optimal cylindrical and full detector veto efficiencies 8 c and as a function 
of the total veto efficiency 8 . Here a = /3 = 0.1 which means that 80% of muons are well 
tracked but 90% are believed to be well tracked. 


The total live time ti in Eq. (5.1) has one very useful feature. It is monotonic in both 
the live times ti in Eq. (3.3) and in Eq. (3.9) of the full detector and cylindrical vetoes, in 
fact it is just a product of a positive powers of the two live times 

ti = [ti (full detector)]^ [ti (cylindrical)]^”^ . (5.3) 

Therefore the tf{ 8 f), Ef{ 8 f), tc{ 8 c), rc( 8 c) and Ec{ 8 c) that minimized the ti of Eqs. (3.3) 
and (3.9) for each individual veto in Sec. 3 also minimize the total ti of Eq. (5.1). This 
greatly simplihes our analysis, as it means that we can continue to use the functions tf{ 8 f), 
Ef{ 8 f), tc{ 8 c), rc{ 8 c) and Ec{ 8 c) that we found assuming perfect tracking. 

Our goal in this section is to obtain 8 f{8) and 8 c{8). When /3 % 0, the former is easily 
obtained from Eq. (5.2) 

Thus we need only substitute 8 f{ 8 ) from Eq. (5.4) into the expression (5.1) for the total live 
time and then choose 8 c{ 8 ) so as to minimize the total live time. 

However the function 8 c{ 8 ) also depends on the parameters a and f3 and so it is not 
easy to display in a hgure. While we have done this optimization for every a and /3, in this 
subsection we display, in Fig. 16, the optimal cylindrical and full detector veto efficiencies 
only for a = /3 = 0.1. This corresponds to the case in which the tracking is successful for 
80% of muons but the experimenter believes that it has been successful for 90% of muons, 
therefore the cylindrical veto is used 90% of the time. Note that in general the optimal full 
detector rejection efficiency is lower than the rejection efficiency for the cylindrical vetoes. 
This is because full detector vetoes lead to a much larger dead time. The maximum possible 
rejection efficiency is 1 — a = 90%, as a fraction a = 10% of the muons are poorly tracked 
unbeknownst to the experimentalist and so the cylindrical veto covers the wrong region. 
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Figure 17: The fractional dead time which is required to obtain a total rejection efficiency 
£ with a = /3 = 0.1 


The corresponding fractional dead time 1 — is plotted in Fig. 17. Although only 80% 
of the muons are successfully tracked, the fact that half of the poorly tracked muons are 
known to be poorly tracked nonetheless allows a rejection efficiency of 85% with a dead time 
of 13%, comparable to KamLAND. 

5.1 The neutrino mass hierarchy 

The optimization of the veto for the neutrino mass hierarchy is now quite simple. Only one 
independent parameter remains to be hxed, the total rejection efficiency £. Similarly to the 
case with perfect tracking, the sensitivity to the hierarchy as a function of £ is obtained by 
restricting Fig. 12 to the pairs {£,ti) on the curve in Fig. 17, corresponding to the maximal 
live time achievable for each rejection efficiency. This yields Fig. 18. Here one can see that 
with a = f3 = 0.1, the optimal rejection efficiency is 78% and it yields a Ax^ of 11.5, as 
compared with 13.8 with no background and 12.4 with perfect tracking. The corresponding 
optimal veto parameters are summarized in Table 3. 

For general tracking efficiency parameters a and /?, one can follow this same strategy 
to optimize the total rejection efficiency £. The results are shown in Fig. 19 together with 
the corresponding fractional live times obtained using Fig. 17. The corresponding statistic 
Ax^ is presented in Fig. 20. This is one of our main results. It shows the dependence of 
the sensitivity to the hierarchy on the tracking. It is evident that a has a greater impact 
upon than f3, reflecting the fact that if it is known that the tracking has failed (as 
occurs for a fraction (3 of muons) then one can nonetheless perform a full detector veto 
which increases the sensitivity to the hierarchy. One can see that the failure to track 10% of 
muons does not greatly reduce the sensitivity to the hierarchy. However, muon bundles are 
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Figure 18: The statistic Ax^ which characterizes the sensitivity to the neutrino mass hier¬ 
archy as a function of the total rejection efficiency S with a = (3 = 0.1 



Figure 19: The optimal veto efficiency and total fractional live time for the determination of 
the mass hierarchy as a function of the tracking efficiency parameters a and (3. The curves, 
from top to bottom, correspond to /? = 0, 0.1, 0.2, 0.3 and 0.4. 


responsible for about 30% of muon events, and so the ability to track bundles does have a 
nonnegligible effect on the sensitivity. KamLAND has several classifications of how reliably 
a given muon was tracked [29]. A tighter reliability threshold corresponds to a larger value 
of (3 and a smaller value of a. The results of our study can be used to optimize this reliability 
threshold. 

We have now optimized the veto strategy, for the determination of the neutrino mass 
hierarchy, as a function of a and (3. For completeness, we will now describe this optimal veto. 
In Fig. 21 one can find the individual veto efficiencies for the cylindrical and full detector 
vetoes. Surprisingly, the individual efficiencies are quite independent of the tracking, one 
always wants a rejection efficiency in the mid to high eighties for the cylindrical veto and 
in the mid sixties for the full detector veto. As these individual efficiencies give the veto 
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Figure 20: The statistic which characterizes the sensitivity to the neutrino mass hier¬ 
archy as a function of the fraction /3 of muons which are known to be poorly tracked and 
the fraction a of muons which appear to be well tracked but in fact are poorly tracked. The 
curves, from top to bottom, correspond to /? = 0, 0.1, 0.2, 0.3 and 0.4. 


parameters, using the results of Sec. 3, this implies that the optimal veto parameters them¬ 
selves are quite insensitive to the tracking abilities of the detector. The optimal parameters 
are shown in Figs. 22 and 23 for the cylindrical and full detector vetoes respectively. 

5.2 Measuring Ou 

The measurement of 612 can be analyzed similarly to the hierarchy determination. The 
precision cT(sin^ ( 26 * 12 )) as a function of S is obtained by restricting Fig. 14 to the pairs 
{S,ti) on the curve in Fig. 17, corresponding to the maximal live time achievable for each 
rejection efficiency. This yields Fig. 24. Here one can see that with a = 13 = 0.1, the optimal 
rejection efficiency is 81% and it yields a (T(sin^ ( 26 * 12 )) of 0.0039 or 0.46%, as compared with 
0.44% with no background or 0.45% with perfect tracking. The corresponding optimal veto 
parameters are summarized in Table 3. 

For general tracking efficiency parameters a and /?, one can follow this same strategy to 
optimize the veto strategy. The corresponding statistic cT(sin^( 26 * 12 )) is presented in Fig. 25. 
One sees that the precision with which 6*12 can be determined is essentially independent of 
the tracking. By performing the same analysis as in the case of the hierarchy, we have also 
found that the optimal veto strategy is also quite insensitive to the tracking, and so for 
brevity we will not repeat the optimal strategy here in general. In the case a = /3 = 0.1 the 
optimal strategy is summarized in Table 3. 
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Figure 21: The optimal cylindrical and full detector veto efficiencies as a function of the 
tracking efficiency parameters a and (3. The curves, from top to bottom, correspond to 
/? = 0, 0.1, 0.2, 0.3 and 0.4. 


Ec (GeV) (meters) (seconds) 



Figure 22: The optimal cylindrical veto parameters as a function of the tracking efficiency 
parameters a and j3. The curves, from top to bottom, correspond to /? = 0, 0.1, 0.2, 0.3 
and 0.4 except for the left panel where the order is inverted 


Goal 

Stat 

S 

ti 

Sc 

Sf 

Ec 

tc 

Tc 

Ef 


MH 

11.5 

78% 

97% 

89% 

66 % 

3.1 GeV 

1.1 S 

2.2 m 

16 GeV 

0.7 s 

sin^(26'i2) 

0.46% 

81% 

95% 

92% 

73% 

2.5 GeV 

1.1 S 

2.3 m 

12 GeV 

0.8 s 


Table 3: The optimal cylindrical veto parameters for the determination of the mass hierarchy 
and 6 i 2 assuming tracking parameters a = /? = 0.1 corresponding to 20% of muons being 
poorly tracked of which half are known to be poorly tracked. The statistics reported are as 
in Table 1. 
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Figure 23: The optimal full detector veto parameters as a function of the tracking efficiency 
parameters a and (3. The curves in the left (right) panel from bottom to top (top to bottom) 
correspond to /? = 0, 0.1, 0.2, 0.3 and 0.4 


cr(sin^(26i2)) (Q'=0.1,y3=0.1) 



Figure 24: The precision with which sin^(26'i2) can be measured as a function of the total 
rejection efficiency £ with a = /? = 0.1 
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Figure 25: The precision with which sin^( 20 i 2 ) can be measured as a function of the fraction 
f3 of muons which are known to be poorly tracked and the fraction a of muons which appear 
to be well tracked but in fact are poorly tracked. The curves, from bottom to top, correspond 
to = 0.1, 0.2, 0.3 and 0.4. For a > 0.25 our maximization did not converge sufficiently 
quickly to be displayed here. 

6 Conclusions 

The next generation of medium baseline neutrino experiments have two main goals, 
the determination of the neutrino mass hierarchy and a precise measurement of 9i2- The 
challenges involved in these two goals are different. The hrst requires a precise measurement 
of the energies of the 15th, 16th and 17th 1 — 3 oscillation peaks in the reactor neutrino 
spectrum [28]. These are small and contain few events, therefore statistical fluctuations 
will be quite important. The second instead requires a measurement of the depth of the 
solar oscillation maximum, a broad feature which extends across the spectrum visible at 
these experiments. This feature includes essentially all of the events in the sample, and so 
it will be measured extremely precisely at these experiments, with an uncertainty which is 
dominated by a systematic uncertainty in the reactor neutrino spectrum [26]. 

As statistical fluctuations play a central role in the determination of the hierarchy but 
not in the measurement of 0 i 2 , the optimal ®Li veto strategies are quite different for these two 
science goals. The hrst is optimized with a very loose veto, which even with perfect tracking 
has an efficiency of only 90%. Of course, given the high rate of muon bundles and showering 
muons [ 6 ] , the tracking will not be perfect and so the efficiency will be appreciably lower. On 
the other hand, as statistical huctuations play a secondary role in the measurement of 6 * 12 , 
a tighter veto is optimal, with an efficiency as high as 93% in the idealized case of perfect 
tracking. This is not a result of a large impact of the spallation background on the precision 
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of the measurement of 612 ■ Indeed, with the 1% uncertainty assumed here in the spallation 
background rate, even with no veto at all the precision in the measurement of sin^( 26 'i 2 ) 
would only increase from 0.44% to 0.47%. On the contrary, it is a result of the fact that a 6 
year run is more than is necessary for a precise measurement of 612 , and so additional dead 
time has little impact on the precision. 

Indeed, a smaller detector with a less PMT coverage could measure 612 nearly as well. 
More to the point, given the central role of systematic errors, we suspect [14] that the most 
precise measurement of 612 at the same price would be given by two smaller detectors, with 
low PMT coverage, at sufficiently distinct baselines so as to break the degeneracy between 
the unknown reactor spectrum and 6 * 12 . Even if this would preclude a determination of the 
hierarchy, it may be the option for RENO 50 which would maximize synergy with JUNO 
and open the door to the measurement of the leptonic CP-violating phase in the years 
beyond [30]. 
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